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Entanglement properties and quantum phases for a fermionic disordered one 
dimensional wire with attractive interactions 
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A fermionic disordered one dimensional wire in the presence of attractive interactions is known 
to have two distinct phases: A localized and a superconducting one depending on the strength of 
interaction and disorder. The localized region may also exhibit a metallic behavior if the system 
size is shorter than the localization length. Here we show that the superconducting phase has 
a distinct distribution of the entanglement entropy and entanglement spectrum distinct from the 
metallic regime. The entanglement entropy distribution is strongly asymmetric with Levy alpha 
stable distribution (compared to the Gaussian metallic distribution), and the entanglement level 
spacing distribution is unitary (compared to orthogonal). Thus, entanglement properties may reveal 
properties which can not be detected by other methods. 
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In the last decade the application of concepts from 
quantum information, such as entanglement entropy 
(EE)[I1, took center stage in understanding physical phe¬ 
nomena in condensed matter physics. One of the reasons 
for this interests is that EE is deeply connected to quan¬ 
tum phase transitions (QPT). The EE quantifies the en¬ 
tanglement in a many-body system by dividing it into two 
regions A and B. For a system in a pure state Ilk), the 
entanglement between regions A and B is measured by 
the EE Sa/b connected to the eigenvalues of the reduced 
density matrix of area A, p a ot H, ps- It is expected 
that non-local properties, such as the EE, may provide a 
different perspective beyond the traditional point-point 
correlations and local order parameters 0. 

Specifically, pA is defined as: pA — TrB|lE')('I'|, where 
the degrees of freedom of region B are traced out. The 
eigenvalues of the matrix are used to calculate the 
EE: 

5^ = -^AflnA^ (1) 

i 

Recently it has been understood that one can utilize the 
information enclosed in the full spectrum of {A)^}. These 
eigenvalues are used to construct the entanglement spec¬ 
trum (ES) {e)^ = — InAf} Q. For one-dimensional (ID) 
systems, the area of the boundary between regions A and 
B is fixed and thus the EE should not depend on the size 
of region A. Nevertheless a logarithmic dependence of the 
form [3-0|' 

= ^In ^sin -f c, (2) 

where La is the length of region A and L is the samples 
length, is expected in the metallic (clean) regime. Several 
years ago Li and Haldane came up with an intrigu¬ 
ing conjecture regarding the connection between the ES 
and the excitation spectrum of a many-body state. They 
suggested that the low-energy ES shows a precise corre¬ 
spondence to the true many-particle spectrum of region 


A. The properties of the EE and the ES may be used to 
identify pha se transitions in disordered many-body prop¬ 
erties 


16 |. 


In this letter we will use the EE and ES in order to in¬ 
vestigate the nature of different phases of fermionic disor¬ 
dered ID systems with attractive interactions. Electron- 
electron interactions in ID systems are parametrized by 


the Luttinger parameter K IT], [18| . For non-interacting 


systems K = 1, while for attractive interactions K > 1. 
When both disorder and interaction are present, an 
extended metallic (with superconducting correlations) 
phase is expected once attractive interactions are strong 
enough, i.e. K > 1.5 [isl - f^ . This stems from the renor¬ 
malization group scaling of the localization length 


^ = (^0)V(3-2iC)^ 


( 3 ) 


where is the non interacting localization length. Thus, 
for K = 1.5 the localization length diverges, and one 
transits from the localized to the extended regime. In¬ 
deed, it has been numerically demonstrated that with 
strong enough attractive interactions in the usual An¬ 
derson model [mnii metal-like behavior emerges, al¬ 
though no evidence of superconducting correlation have 
been numerically demonstrated. 

For disordered system one must consider the EE be¬ 
havior over an ensemble of different realizations of dis¬ 
order. It has been demonstrated in Ref. ll| that the 


median EE for La < ^ follows the metallic logarithmic 
behavior (Eq. ([2])), while for La > ■C it saturates. In 
principal this could be used to decide in what phase (lo¬ 
calized or metallic) the system is in [l^. Nevertheless, in 
a realistic numerical study this strategy is fraught with 
problems, since the localization length grows fast as func¬ 
tion of K (Eq. ([2])), and easily outgrows any finite sys¬ 
tem length L. Once ^ ^ L a finite system will show 
metallic behavior although it is in the localized regime. 
For brevity, we shall refer to the K > 1.5 regime as su¬ 
perconducting, and to the finite sample K < 1.5 regime 









as localized or metallic according to whether L > ^ or 

i>L. 

In this Letter we will show that the full distribution 
of the EE shows a distinct behavior between the metallic 
and superconducting regime, although the median EE is 
essentially identical in both regimes. The EE distribu¬ 
tion changes from a Gaussian in the metallic regime to 
a very asymmetric Levy alpha stable distribution with 
“fat tails” in the superconducting regime. The ES level 
spacing (ESLS) distribution fits the Gaussian orthogo¬ 
nal distribution (GO E) exp ected from interacting many 
particle system [3, [25l - [30| in the metallic regime, while 

it changes to a Gaussian unitary distribution (GU E) as¬ 
sociated with superconducting excitations [8lL l82j| in the 
superconducting regime. Thus, the EE and ESLS distri¬ 
butions are able to characterize the phase of the system, 
where other methods fail. 

In this letter we consider a spinless ID electrons wire of 
size L with attractive nearest neighbor interactions and 
on-site disordered potential. The system’s Hamiltonian 
is given by: 

L L-l 

^ ( 4 ) 

1=1 i=i 

1 I 

+ ^ “ 2 ^’ 
i=i 


where ej is the on-site energy, which is drawn from a uni¬ 
form distribution [—IT/2, IT/2], c] is the creation opera¬ 
tor of a spinless electron at site j , and < = 1 is the hopping 
matrix element. The interaction strength is C/ < 0, and 
a background is included. For the non-interacting An¬ 
derson model the system is localized with a localization 
length ^0 ~ 105/IT^ [s^. Here the Luttinger parameter 
K{U) = 7r/[2cos“^(—17/2)] For non-interacting 

electrons K(U = 0) = 1. For attractive interactions 
K > 1 and ^ increases as U becomes more negative. For 
U = —1, K = 1.5 and the localization length according 
to Eq. diverges. Thus, below 17 < —1 the system is 
expected to be delocalized. At 17 = — 2 it goes through 
another phase transition to a phase separated state and is 
insulating again. Indeed numerically 3) 22 - 24 j, this sys¬ 
tem is known to exhibit extended behavior for a range of 
attractive interaction strength centered around U = —1.5 
and not too strong disorder IT < 1.5 

The density matrix renormalization group (DMRG) 
s a very accurate numerical method for calculat¬ 


ing the ground state of disordered interacting ID system 
and for the calculation of the reduced density matrix. 
We calculate the EE for three length L — 800, 700,1100 
and different values of La = 10,20,... ,L — 10, for 
400, 200,100 realizations of disorder for the correspond¬ 
ing system length. Specifically, we calculate the normal¬ 
ized EE of the j-th realization at a given La, Sj^La) = 
Sj{LA)/{S{L a)), where {S{La)) is the average EE over 



FIG. 1: (Color online) The distribution P{s) of the normal¬ 
ized EE (see text) for different length L = 300 (black, circles), 
L = 700 (red, squares), L = 1100 (green, diamonds), and dis¬ 
order strength W — 0.3 (empty symbols), IT = 0.7 (full sym¬ 
bols). A fit to a Gaussian with a width which depends on the 
disorder, g{W), is depicted by the continuous lines, (inset) 
The median EE as function of La. The symbols correspond 
to the numerical results, the curves to Eq. ©• 


the different realizations. Since the distribution of the 
EE is very similar for different values of La as long as 
La is not to close to the edge we accumulate the dis¬ 
tribution of the normalized EE, P{s), in the range of 
L/4 <La< 3L/4. 

Let us begin by discussing the EE distribution for 
U = —0.7 {K{U = —0.7) = 1.8) for which the system 
is in the metallic regime, i.e., localization length much 
larger than sample length. In Fig. [I]we present the dis¬ 
tribution for two values of disorder (IT = 0.3, IT = 0.7), 
corresponding to the non-interacting localization length 
Co(IT = 0.3) - 1200 and ^o(VE = 0.7) - 200, thus 
^(IT = 0.3,1/ = -0.7) - 5 • 10^ and ^(IT = 0.7, U = 
—0.7) ^ 6 • 10®. First, we examine the median value of 
the EE as function of La (here we could use the average 
which is almost equal to the median, but for the sake of 
uniformity with the upcoming results we use here the me¬ 
dian) and compare it with the expression of the EE for a 
clean system described in Eq. As can be seen in the 
inset of Fig. [l]it fits quite well. Thus, the median EE fol¬ 
lows closely the expected behavior of a clean (metallic) 
system. The distribution P{s) for all three length and 
two disorder strength is plotted in Fig. [TJ As has been 
discussed in Ref. [38|, for L <C we expect the distribu¬ 
tion to be Gaussian (i.e. exp(—(x — l)^/2cr^) 

centered at the average. This is indeed seen in Fig. [I](up 
to a slight skewness of the tails), as well as the fact that 
the width, a of the Gaussian is almost independent of L. 
On the other hand it is clear that the disorder strength 
IT does affect a. As might be expected, The weaker is 
the disorder, the narrower is the width of the Gaussian. 

The dependence of the Gaussian width on IT is ex¬ 
amined in Fig. 121 We keep the length and interaction 
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FIG. 2: (Color online) The distribution P{s) of the normal¬ 
ized EE for a fixed wire length (L = 300) and different disor¬ 
der strength: W = 0.3 (black circles), W = 0.5 (red squares), 
W = 0.7 (green diamonds), W = 0.9 (blue triangles). A fit 
to a Gaussian of width, a{W) (depicted in the inset) is indi¬ 
cated by the continuous lines, (inset) The Gaussian width a 
as function of the disordered strength square W^. 

strength fixed (L = 300, U = —0.7) while varying W. 
Even for the strongest disorder ^{W = 0.9, U = —0.7) ^ 
2 • 10®, is much larger than the sample length. P{s) re¬ 
mains Gaussian for all disorder strength. Moreover, as 
can be seen from the inset a{W) <x. . 

What happens to the EE distribution in the supercon¬ 
ducting regime? Specifically, we concentrate on the ex¬ 
tended regime with U = —1.5 {K{U = —1.5) = 2.3), 
different values of disorder W = 0.3,0.5,0.7,0.9 cor¬ 
responding to a non-interacting mean free path 
1200,400,200,130 and different sample length L = 
300,700,1100. This regime of the parameter space is 
deep in the superconducting regime. We present the dis¬ 
tribution P(s) for all length and interaction strength in 
Fig. [21 It is obvious that P(s) is completely different 
than in the metallic regime (Figs. I1I2I) . In all cases the 
distribution is strongly asymmetric. The typical value of 
the EE is close to its maximum value, and the probabil¬ 
ity of measuring an EE larger than the typical value is 
rather small and falls off abruptly. On the other hand, 
there is a high probability for measuring values of the 
EE below the typical value. For smaller values of the EE 
the distribution has a very long tail. Another difference 
is the strong dependence on the sample length L even at 
the same value of disorder W. 

First, one concludes that similarly to the case of the EE 
of a localized system (L > ^) [ll|, [ij, , the average EE 
is not a correct description of the superconducting regime 
typical EE. The EE is more appropriately represented by 
the median value. In the inset of Fig. [3] the median value 
of the EE as function of La is compare to Eq. ([2]). As 
can be seen it fits quite well. 

As can be seen in Fig. [U the distribution is uni¬ 
versal and may be scaled by the function P{u), here 



FIG. 3: (Golor online) The distribution P{s) of the normal¬ 
ized EE deep in the superconducting regime {U = —1.5) 
for different strength of disorder (W = 0.3, black symbols, 
W = 0.5, red symbols, W = 0.7, green symbols, W = 0.9, 
blue symbols), and sample length L = 300 (circles), L — 700 
(squares), L = 1100 (diamonds), (inset) The median EE as 
function of La- The symbols correspond to the numerical 
results, the curves to Eq. ©. 

M = l-|-(s—1)* Pm&x, and the distribution is normalized 
P(u) = P(u)/Pmax, where Pmax IS the value of P(s) at 
the maximum. In the inset of Fig. [?] , the values of Pmax 
vs. ^o/L is depicted. Up to the ballistic regime, (i.e.,for 
^0 P), a linear relation seems to hold. The distribution 

P{u) may be described rather well by a Levy alpha stable 
distribution, which is a natural extension for the central 
limit theorem to the case that the identically-distributed 
random variables have no finite variance, and is used to 
describe distributions with “fat tails”. The Levy distri¬ 
bution 

OO 

f{x,a,l3,j,6) = ;^Re j (5) 

0 

(where $ = tan( 7 rQ!/ 2 ), except for a = I where $ = 
— ( 2 / 7 r) log(t)), is defined by four parameters. The sta¬ 
bility index 0 < a < 2 characterizes the asymptotic be¬ 
havior of the tails (except for a = 2 ), the skew¬ 

ness parameter —I > /3 < 1, 7 is a scale factor and 
S controls the location of the maximum [^ . In gen¬ 
eral / is not analytical, except for special cases, such as 
fix, a = 2 ,/ 3,7 = (tI\/ 2,5 = x), equal to a Gaussian of 
width a. We plot /(u, a = 1, (3 = —1 ,7 = 0.285, S = 0.8), 
and fiu,a = 1.2,/3 = - 1,7 = 0.28,J = 0.52) in Fig. 
m The skewness is clearly maximal (/3 = —I), while 
fitting a depends very much on the tail region, which 
has the largest numerical uncertainty. Nevertheless, the 
main part of the distribution is evidently fitted well by 
1 < a < 1.2. 

Thus, although the median of the EE (or its average) 
does not give us a clear signature whether the system 
is in the metallic regime or in the superconducting one. 
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FIG. 4: (Color online) The data presented in Fig. [3]scaled by 
the parameter u (see text). The continuous (dashed) black 
curve is the Levy stable distribution with a = l,/3 = — 1,7 = 
0.285,5 = 0.8, (a = 1.2, P = - 1,7 = 0.28,5 = 0.52). (inset) 
The maximum of the distribution P{s) (Pmax) depicted in 
Fig. [3] as function of S,o/L. A linear behavior in the regime 
^o/L < 2 is indicated by the line. 


the distribution can differentiate between the regimes. 
Further work is needed in order to follow how the distri¬ 
bution changes as one goes through the phase transition 
and is the transition point characterized by a special dis¬ 
tribution. A hint may be found in Ref. [iSj, where the 
variance of the EE is plotted through the transition. It 
seems that the variance becomes smaller in the metal¬ 
lic side as one approaches the transition and then grows 
again in the superconducting regime. This is in line with 
a crossover from a Gaussian distribution to a fat tail one. 

In calculating the ESLS one must take into account 
that different sectors of the reduced density matrix with 
a fixed number of electrons do not couple with each other. 
Therefore, for each eigenstate of pA one should calculate 
both Xi and (the number of particles in the region 
A). Calculating N:^ does not add to the complexity of 
the DMRG code 0. Thus, each eigenvalue of the reduced 
density matrix has two indexes X ^'^, which are translated 
into the entanglement Hamiltonian eigenvalues Q = 
— ln(Af^^). It makes sense to calculate the ESLS only be¬ 
tween eigenvalues belonging to the same number sector. 
Thus, and omegaf^ = Af^/(Af'^). 

The average spacing for U = —0.7 and U = —1.5 for 
Na = L/2 = 350 as function of i is depicted in the inset 
of Fig. [5] (Aj' ) only slightly depends on U. There 
are several anomalously large values appearing for the 
first few spacings. The distribution P{uj) is accumu¬ 
lated over all low-lying spacings < 20), which are 

not anomalously large (i.e., (A^*^) < 1) for all values of 
Na, and L/4 < La < 3L/4. The distributions for the 
metallic {U = —0.7) and superconducting {U = —1.5) 
regimes are presented in Fig. [S] It can be seen that 
for the metallic case the numerical results fit quite well 



FIG. 5: (Color online) The ESLS distribntion P{uj) as func¬ 
tion of the normalized level spacing oj. The numerical results 
for the metallic {U = —0.7, black circles) and superconduct¬ 
ing ([/ = —1.5, red squares) regimes for L = 700 are indi¬ 
cated. The black (red dashed) line correspond to the GOE 
(GUE) distribution. Inset: The average level spacing {A^^} 
as function of level number for both interactions. 


the GOE distribution (Pgoe(s) = (7rs/2) exp(—7rs^/4)) 
expected for the excitations of interacting many particle 
systems @,i,[2ii3- On the other hand, in the super¬ 
conducting regime the distribution follows the GUE dis¬ 
tribution (F’gue(s) = (32s^/7r^) exp(—4s^/7r)). Indeed, 
superconducting correlations lead to exotic statistics of 
the excitations close to the Fermi energy [^, which are 
quite similar to GUE for higher energies. Thus, also the 
ESLS distribution shows a significant difference between 
the metallic and superconducting phases, although for 
the length considered here (L = 700 <C ^ = 6 • 10^) both 
regions are extended. 

To conclude, the EE and ESLS distributions show 
a distinct behavior, depending if the system is in a 
metallic (system length much smaller than the localiza¬ 
tion length) or a superconducting (localization length di¬ 
verges) regime, for a disordered one dimensional spin¬ 
less electron system with attractive interactions. While 
the metallic regime shows a rather expected behavior, 
e.g., Gaussian distribution of the EE and GOE ESLS 
distribution, the superconducting EE shows an asymmet¬ 
ric distribution and GUE ESLS. Thus, the entanglement 
properties encode details of the underlying phase of the 
system which may elude other measures, which may be 
useful for detecting new phases for different systems. The 
emergence of a “fat tail” distribution for the EE in the 
superconducting regime is fascinating and deserves fur¬ 
ther study. 
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